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We improve previous calculations of the CMB spectral distortions due to the decay of primordial 
magnetic fields. We focus our studies on causally generated magnetic fields at the electroweak 
and QCD phase transitions. We also consider the decay of helical magnetic fields. We show that 
the decay of non-helical magnetic fields generated at either the electroweak or QCD scale produce 
p and y- type distortions below 1CU 8 which are probably not detectable by a future PIXIE-like 
experiment. We show that magnetic fields generated at the electroweak scale must have a helicity 
fraction /* > 10~ 4 in order to produce detectable p- type distortions. Hence a positive detection 
coming from the decay of magnetic fields would rule out non-helical primordial magnetic fields and 
provide a lower bound on the magnetic helicity. 


I. INTRODUCTION 

Magnetic fields are observed throughout the cosmos; 
from galaxies at high and low redshifts m, in galaxy 
(super)clusters [UG>], an d in the voids of the large scale 
structure [5] • Although still under debate, the difficulties 
of astrophysical mechanisms in explaining such observa¬ 
tions promote the idea of a primordial origin for magnetic 
fields i.e. fields generated in the early Universe before 
structure formation. There are a number of theoretical 
mechanisms proposed to generate primordial magnetic 
fields, such mechanisms for example involve inflation [7], 
first-order phase transitions [8], or vorticity generation 

nmm. 

If primordial magnetic fields were indeed generated, 
they would suffer decay on small scales due to magneto¬ 
hydrodynamics (MHD) effects [TIMT5j . The dissipation of 
magnetic fields injects energy into the plasma, if this oc¬ 
curs in the early Universe when the CMB is being formed 
distortions to its black-body spectrum can be gen¬ 
erated [T5U25] . Distortions to the CMB black-body spec¬ 
trum come in different types depending on the epoch of 
energy injection. At relatively early times a p-type dis¬ 
tortion can be generated, where the photon Bose-Einstein 
distribution develops a non-vanishing chemical potential. 
Whereas at later times a Compton y-parameter can be 
generated giving a y-type distortion. Mixed distortions 
are also possible [25]. Any mechanism which injects en¬ 
ergy into the plasma in the early Universe has the po¬ 
tential to generate such distortions. Other mechanisms 
in the early Universe include the dissipation of primor¬ 
dial acoustic waves 01 HE], decaying or annihilating relic 
particles (2TH25] , and the evaporation of primordial black 
holes or cosmic strings [30, 31j. 

In this paper we consider magnetic fields generated by 
some causal process in the early Universe. In particu¬ 
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lar we consider the possibility of magnetic fields gener¬ 
ated by first-order phase transitions at either the elec¬ 
troweak (EW) or the quantum chromodynamics (QCD) 
scale [8]. The spectrum on large scales for magnetic fields 
generated in such a process is highly constrained due to 
causality reasons [32]. The slope of the large scale spec¬ 
trum also determines the subsequent MHD evolution of 
magnetic energy and coherence scale [T211151155] . As the 
small scales dissipate into heat in a turbulent plasma, the 
peak of the spectrum moves down along the large scale 
spectrum [H [H [51 [55]. The amplitude of magnetic 
helicity also determines the evolution of magnetic fields. 
Helicity conservation in the early Universe slows down 
the decay rate of fully helical fields compared to non- 
helical fields and leads to an inverse cascade of energy 
from small scales to large scales mm- Hence causally 
generated magnetic fields and their initial helicity frac¬ 
tion can be constrained by CMB spectral distortions. 

In Section |TT] we use the results of Refs. [nun [55] to 
determine the decay rates of magnetic fields as functions 
of their initial strength, coherence length and initial he- 
licity fraction. From these decay rates, in Section uni 
we calculate the CMB spectral distortions generated and 
analyse the parameter space. We conclude in the final 
Section. 


II. THE DECAY OF MAGNETIC FIELDS IN 
THE EARLY UNIVERSE 


In order to calculate the CMB spectral distortions we 
need to know the energy injected into the primordial 
plasma. For decaying magnetic fields, the energy injec¬ 
tion rate is given by [20] EU 


dQ _ __ _ 4 d p B 

d t d t 


(1) 


where ps is the comoving energy density i.e. pB = aPps- 
The averaged magnetic energy density is obtained by in¬ 
tegrating over the local energy density ub = B 2 /87r, here 
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B = a 2 B is the comoving magnetic field, 

Pb = y Ju B dr=-^~ J |B(k)| 2 dk = pj M k dk, (2) 


where p is the total comoving energy density and Mk is 
the magnetic spectral energy. Here we emphasise that the 
spectrum on large scales, parametrized by M*, oc k n+2 : 
for causally generated magnetic fields is constrained by 
n > 2 [32]. The most shallow, and expected, scaling 
n = 2 is also confirmed by numerical simulations giving 
Mk oc k 4 [531155] . Assuming that the magnetic energy is 
concentrated at the Integral scale (index /), which defines 
the peak of the spectrum in Fourier space, we can write 
Pb = P J kMfrd In k ~ pAqMj, adopting the conventions 
of Ref. [53] where the wave vector k is also comoving. 
Since the photon energy density scales as p 7 oc a~ 4 , we 
can write 


_ 1 dQ __ _i dps _ 1 d{pk I M I ) 
Pi dz - Pl,0 dz - P 7,0 dz 


(3) 


where ao = a(To) = 1 today. 

The magnetic field strength Bj and coherence length 
Li (identified as 2n/ki) evolve during the radiation dom¬ 
inated era due to turbulent MHD effects [T3}[5TH551. In 
the following we quote the results from the detailed an¬ 
alytical and numerical analysis in Refs. PH[35][55]. The 
evolution of the field strength and coherence length de¬ 
pends on the state of the plasma and on the properties of 
the magnetic field, in particular its hclicity. The average 
helicity density is given by 


hs = 


1 

V 


(A ■ B)dr = p 


Bkdk , 
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where B = V x A. For a magnetic field with helicity, on 
any given scale k there is a realizability condition given 
by \Bk\ < 87 rMk/k, where the helical spectrum Hk is de¬ 
fined following the conventions of Ref. [55] , From the 
above we can define / = k'Hk/8nMk as the helicity frac¬ 
tion, where / = 0 for the non-helical case and / = 1 
for the maximally helical case. The helicity density is a 
useful quantity since it is conserved in the early Universe 
h B — const, when the conductivity er = l/Anr] —► 00 [30] . 
and the conservation of magnetic helicity determines the 
evolution of the magnetic field strength and coherence 
length. 

In the turbulent regime, where kinetic Reynolds num¬ 
bers are large R e ^> 1 , the general decay law for the 
magnetic energy is Mj oc a - 2 ( n + 2 ) /O+s) anc j f or the co¬ 
moving integral scale Lj <x a 2 /( n + 5 ) [TH 3U [35], hence 
Bj oc a~ 5 / 7 , ki oc a -2 / 7 and p B oc a -10 / 7 , where n = 2 
is used due to causality constraints for the large scale 
spectrum PU EH EH Eg. These decay laws for non¬ 
helical magnetic fields, where the helicity fraction / <1 
if not zero, are obtained through analytical considera¬ 
tions in Refs. nans os] and confirmed numerically in 
Ref. P2]- If the magnetic field has non-zero helicity, the 
helicity density in this case would grow as Hf oc a 2//7 , 


and eventually the magnetic field would become fully he¬ 
lical / = 1. The conservation of magnetic helicity en¬ 
sures the relation p B oc B 2 oc kj. The decay of magnetic 
fields slows down in the fully helical case and the de¬ 
cay rate becomes independent of the shape of the large 
scale spectrum n. For the maximally helical case we find 
Bj oc a -1 / 3 , ki oc a -2 / 3 and p B oc a -2 / 3 . We can sum¬ 
marise the above turbulent damping (TD) decay laws by 

_np+5 2 

(TD) : p B oc a *>+ 7 , Li oc a?+ T , (5) 

where p = 0 for non-helical fields and p = —4 for maxi¬ 
mally helical fields. 

Here we note that an apparent inverse transfer of mag¬ 
netic energy has been numerically observed for the non¬ 
helical case. This effect leads to a weaker evolution 
for non-helical magnetic fields in the turbulent regime 
Bi oc a -1 / 2 , ki oc aT 1 / 2 and p B oc a -1 , i.e. p = —3 
above, giving Bj oc Lj 1 [33] (see also Ref. 32])- This 
is potentially a very interesting and exciting new devel¬ 
opment in turbulent MHD. However this is a numerically 
observed effect under the conditions of high resolutions, 
and magnetically dominant turbulence [361 . The condi¬ 
tion of magnetically dominant turbulence is perhaps not 
satisfied in the early Universe. Magnetogenesis at first or¬ 
der phase transitions typically produce a lot of turbulent 
kinetic energy. The generated magnetic field, through 
dynamo action, comes into equipartition with the kinetic 
energy, but is unlikely to dominate over the kinetic en¬ 
ergy, see for example Ref. 0. Furthermore, in the study 
by Ref. (56] it seems that the inverse transfer is less effi¬ 
cient for large Prandtl numbers, but the Prandtl numbers 
in the early Universe are huge. In any case we will also in¬ 
vestigate the CMB spectral distortions considering these 
scalings to complete the study. 

In the viscous regime, where R e < 1, during particle 
diffusion, the rapidly growing particle mean free path 
means that the dissipative time scale is increasing faster 
than the Hubble time if -1 . This prevents further dissi¬ 
pation of magnetic energy, and so the magnetic energy 
freezes-out and remains constant fi~2 

(VF) : Pb — const., L/ ~ const., (6) 

The above evolution in the viscous freezing (VF) stage 
occurs regardless of magnetic helicity. However, as the 
mean free path increases further particles begin to free- 
stream out of overdensities. In this case the interaction 
between the fluid and the background is described by a 
drag force, which is a decreasing function of time. This 
leads to a situation where the magnetic dissipative time 
scale can become smaller than the Hubble time once 
again and magnetic energy decay can start again ram. 
The drag force due to the free-streaming of neutrinos be¬ 
comes ineffective at approximately the freeze-out of neu¬ 
trinos, i.e. when their mean free path becomes larger 
than the Hubble scale. The dissipation due to photon 
drag is efficient until photon freeze-out at T ~ 0.26 eV 
when recombination commences. The decay of mag¬ 
netic fields due to viscous photon free-streaming is given 
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by jfcj oc a ~ 3 /( n + 5 ) and p B cx a -3(n+3)/(n+5) [jg. With 
n = 2 we find I?/ oc a -15//14 , ki oc a -3 / 7 and ps oc a -15 / 7 , 
which is a faster decay than in the turbulent case, and 
faster than in the case of a maximally helical field. 
For maximally helical fields the relation p B oc kj again 
holds, and the decay due to free-streaming photons is 
Bj oc a -1 / 2 , ki oc a -1 and p B oc a -1 . We can summarise 
the above viscous damping (VD) decay laws by 

. . _o pd~5 3 

(VD) : p B oc a p+ 7 , T/ oc a* > + 7 , (7) 


for velocity fluctuations v™ s correlated on some physical 
scale l. Here p* and g 1 are the total and component 
number of effective relativistic degrees of freedom. We 
apply the above from the time of neutrino decoupling at 
T^dec ~ 2.6 MeV to the time of matter-radiation equality 
at T eq — 1 eV. In this epoch, the comoving photon mean 
free path is given by @3] 


77 

^mfpjC 


VTyJnl^ + nl ’ 
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where p = 0 for non-helical fields and p = —4 for maxi¬ 
mally helical fields. 

In order to calculate the CMB spectral distortions due 
to decaying magnetic fields, we have to know in which of 
the above phases of evolution the p and y eras correspond 
to. For this we must understand when these phases begin 
and end, which is what we do in the next Section. 


A. End of the turbulent damping stage 

When the magnetic fields are first generated, for ex¬ 
ample by a first-order phase transition, the magnetic en¬ 
ergy is expected to come into equipartition with the ki¬ 
netic energy [see Ref. [SS] f° r a detailed analysis of this 
mechanism]. At the time of the EW or QCD scale it 
can be shown that the kinetic Reynolds numbers are 
very large in this case and the plasma is highly turbu¬ 
lent. Hence magnetogenesis ends with the plasma in a 
turbulent state. We now calculate Xe 0 t, the temper¬ 
ature corresponding to the end-of-turbulence. Turbu¬ 
lence ends when the kinetic Reynolds number decreases 
to R e ~ 1, hence the end-of-turbulence temperature is 
obtained through 


Re (Teot) ~ 1 • (8) 


where ut = 8ira 2 /3nil is the Thomson cross section, 
a ~ 1/137 is the fine structure constant and m e is the 
electron mass. The number densities n pair and n e of er- 
pairs and free electrons respectively are given by [431] 


n 


pair 


( 2 m e T \ 
(—) 



Tl e — X e 


flbPo 

TOpr 




( 11 ) 

( 12 ) 


where m pr is the proton mass, the baryon fraction and 
present day density product is flbPo — 1.81 x 10 -12 eV 4 
[44] . T 0 — 2.725 K is the present day photon tempera¬ 
ture and the ionization fraction is X e = 1 in the radiation 
dominated era. Following the usual assumption that in 
the turbulent regime R e (Lj) 1 there is equipartition 
between magnetic and kinetic energy on all scales up to 
the integral scale, i.e. 


/ 6rms\2 

(u™ 8 ) 2 = k ’ for k>ki. (13) 
47T/9 

Numerical simulations show that there is almost exact 
equipartition i.e. T ~ 1 for non-helical fields, which 
is slightly reduced to T « 10 _1 for maximally helical 
fields [12]. Therefore we find on the integral scale 


At the time of the EW or QCD scale the plasma viscos¬ 
ity is generated by neutrinos, since at this time they are 
the particles which are the most efficient at transporting 
momentum and heat m- But as the Universe expands 
and cools the neutrino mean free path increases which in¬ 
creases the plasma viscosity. Therefore, due to neutrinos, 
the plasma goes from a turbulent state to a viscous state 
up until the neutrinos decouple at T^ec — 2.6 MeV. The 
evolution of magnetic energy and its coherence length in 
the turbulent stage due to neutrinos followed by the vis¬ 
cous stage before neutrino decoupling can be well approx¬ 
imated by only considering a turbulent damping stage 
throughout that epoch |T2] . In our calculations we use 
this simplifying approximation. 

After neutrino decoupling we are in the photon era, 
when photons generate the plasma shear viscosity r) s . In 
this case the Reynolds number is given by m 


Re(l) 


s 


5 ff *(T) v\™l c 
g 7 Gp.cCn’ 


(9) 


Re(Li, T ) ~ Vf ^ L/ ( 2 , (14) 

57 mfp,c V P J 

which is to be evaluated from the time of neutrino 
decoupling. To connect the magnetic energy density 
with the coherence length we use the approximation dis¬ 
cussed below Eq. (|8|. Hence we can use Eq. ([5]) to find 
Li\Tpb = Ri,*\/pB,* (T/T*) (p+3)/(p+7) , where the index 
**’ denotes the epoch of magnetogenesis. With the above 
the Reynolds number is 


R e (Lj,T) ~ V2IW/3 


5 g* A 


max 

B,* 


9v,l ^mfp.c 


yV 


(15) 


where we have defined e = p Bf */p lt o — p B ,*/p- Here, 
the maximal magnetic energy e = 1 corresponds to 
its = p/2, i.e. a maximum magnetic field strength of 


BZT = \f^~P — 3 x 1(T 6 G , 


V, 


(16) 
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where the radiation here is taken to be the CMB 
photons 12] • In the above we have also defined 
P = Ab,*/A^ x , and we have identified the integral scale 
with the magnetic field coherence length A b,*- For 
magnetic fields generated at a time during the radiation 
dominated era (in contrast to inflationary magnetogene¬ 
sis), the basic constraint on the coherence length is the 
horizon size at the time of magnetogenesis 


where pb is the baryon density. The Alfven velocity is 
given by v\ ~ 2ps/Pb, since at this time the photons are 
decoupled from the fluid. Magnetic dissipation due to 
photon drag is shown to be efficient until photon decou¬ 
pling |T?| . 

Hence we can solve Eq. (19) for T V( i, the temperature 
at the start of the viscous damping due to free-streaming 
photons. We find 


A b,* < A 


max 
B ,* 


l 

aH * 


(17) 


The horizon size is 2 x 10 -10 Mpc and 3 x 10“' Mpc at 
the electroweak and QCD phase transitions respectively. 
To estimate the integral scale at the time of magnetogen¬ 
esis T/,*, we can assume that turbulence is effective such 
that the Alfven eddy-turnover time Ia,* is equal to the 
Hubble time [T2] he. 


i a ,* — 


L 


,rms 

; A,L 


i 

a.H 


(18) 


where v ™£ = B™ s / y/3itp. With the above we find that 
A b,* = Lj : * = y/e/aH\ *, hence we can set p = yfe in 
Eq. ( fib] ) above. 

The end-of-turbulence temperature, obtained through 
R e (Teot) ~ 1, can only be determined numerically; we 
cannot analytically invert the function in Eq. (15) due 


to the exponential in the photon mean free path. We 
find that T Eo t — 2 x 10 4 eV and 2 x 10 3 eV for non¬ 
helical magnetic fields generated at the EW and QCD 
phase transition, T* ~ 100 GeV and 200 MeV respec¬ 
tively, with p = 1 and e = 1. Hence the //-era is within 
the viscous regime, since the //-era commences at around 
~ 470 eV. However, for maximally helical fields, it 
is possible that the plasma is still in a turbulent stage 
within the //-era. The point here is that the magnetic 
field decay could be in the turbulent, viscous freezing or 
viscous damping regime depending on the initial condi¬ 
tions. This can be seen in the evolution plots in Figure]!] 


T 3 9 / 90 m P m pl m 2 e T 0 ( T EoT \ 

vd ~ i6^ 2 a 2 y 3 *, vd x eM n bPo 1 ’ 

which is valid for T < 10 4 eV when the photon mean free 
path can be well approximated by /[[ ifp ~ l/arfie- 

The evolution of magnetic energy due to MHD turbu¬ 
lence occurs approximately until the time of recombina¬ 
tion (index ‘rec’), when the field configuration falls on 
the line given by mm 

B Arec ~8xlO- 8 ^f^G. (22) 

’ Mpc v ’ 

This line corresponds to the largest eddies being pro¬ 
cessed at recombination l/(aff)| rec ~ X/va with va the 
Alfven speed B3 05]. Beyond this epoch, the evolu¬ 
tion of the field strength and coherence length essen¬ 
tially ceases, with only a logarithmic scaling [12] . and the 
magnetic fields become frozen into the plasma. Hence, 
fields generated during the radiation era will evolve to 
fall on the above line at recombination, which are also 
the values that will be observed today B\ rec ss Bq and 
As,rec ~ A b, since the field strength and coherence length 
do not evolve significantly in the matter dominated Uni¬ 
verse. 

We have now all the necessary ingredients to charac¬ 
terise the full evolution of magnetic fields in the radiation 
dominated epoch. Exemplary evolution histories of mag¬ 
netic energy for varying initial conditions can be seen in 
the plots of Figure [l] where different initial conditions 
A b,*, and /* lead to different histories, and hence 
different decay rates during the // and y-eras. 


B. Start of the viscous damping stage 


The magnetic energy and coherence length stop evolv¬ 
ing in the viscous regime when particles are diffusing 
lmf P ,c ^ L [see above Eq. §]■ The viscous damp¬ 
ing stage occurs when photons begin to free-stream 
^mfp,c L. The start of the viscous damping stage (in¬ 
dex ‘vd’) is determined by the condition [12] 

Tvisc.free (T vd ) = H- 1 (T vd ) , (19) 


where the viscous free-streaming time-scale is given by 
Arise.free = a 7 L 2 /v\. Here, the drag term on the fluid 
due to the occasional scattering of photons with the fluid 
particles is given by [12] 


4 1 Pa 
3 ^mfp Pb 


( 20 ) 


III. SPECTRAL DISTORTIONS FROM 
DECAYING MAGNETIC FIELDS 

At high temperatures, corresponding to :>2x 10 6 , 
the black-body spectrum of the CMB [12] is formed from 
bremsstrahlung and Double-Compton scattering [see for 
example Ref. [20] and references therein]. As the red- 
shift drops below 2 x 10 6 these interactions become in¬ 
efficient at restoring the black-body spectrum if addi¬ 
tional energy is injected into the plasma and distortions 
could be imprinted from then on. In the early stage 
2 x 10 6 > z > 5 x 10 4 the elastic-Compton scattering is 
efficient enough and the spectral distortion comes in the 
form of a non-vanishing chemical potential //. This //- 
type CMB spectral distortion is generated if thermal en¬ 
ergy is injected into the plasma during the //-era defined 
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and Yp = 0.24 is the primordial helium mass abundance 
[46] . The solution to Eq. (23) is given by [20] 




dz d Q 
— a“ ex P 

p 1 d^ 



(25) 


where 

1.97 x10 s / n b h 2 
ZDC = (a _ IJi_\f V0-0224 ) 

V 1 2 0 . 24 / 


(26) 


and where z; = 2x 10® and z en d = 5 x 10 4 define the 
start and end of the /z-era. 


A. CMB /z-type distortions from decaying 
magnetic fields 

From the full evolution history of magnetic fields shown 
in Figure[l]we can see that, in most cases, the plasma is in 
the viscous regime during the /z-era. This is true in most 
cases except for fields generated at the QCD scale if the 
initial helicity fraction is large enough. For analytical 
purposes let us first consider the viscous damping law, 
from which we can write 


PB\~) - pb 




1 + 


o P+5 
6 P+7 


(27) 


FIG. 1. In plots (a) and (b) we show the evolution of mag¬ 
netic energy from the time of magnetogenesis (at the EW 
and QCD scales respectively) to recombination. In plots (a) 
and (b) the (solid, blue) lines from top to bottom correspond 
to initial helicity fractions /* = {lO -3 ,10~ 4 ,10~ 6 , < 10~ 14 } 
and /* = {10 —1 , 10 —4 ,10 -6 , < 10~ 14 } respectively. The max¬ 
imally helical case /* = 1 (solid, red) is also shown. We 
also show (dashed, gray lines) the evolution of non-helical 
magnetic fields with an inverse transfer of energy (see dis¬ 
cussion below Eq. §). In all the plots shown we con¬ 
sider e = pB,*/p~t,o ~ 1, which corresponds to an initial field 
strength of Bx,* — 3 x 10~ 6 G, with e < 1 the evolution his¬ 
tory also changes. 


above. The rate of change of the chemical potential p is 
determined by (TfiUIfi] 


dp p 1.4 ^dQ 

df fDc(~) 3 d t 1 


(23) 


where dQ/dt is the energy injection rate. Here, the time 
scale for Double-Compton scattering is 


where pp,v d indicates the magnetic energy at the start of 
the viscous damping stage. Here we also use a oc 1/T for 
the photon temperature and T = To (1 + z). Since the 
magnetic energy is frozen-out from the time of the end- 
of-turbulence (EoT) to the start of the viscous damping 
stage (vd) we can set ps,v d — Pb, EoT, hence we find 


-id Q 

P7 dz 


_i d pb _ 

" P7 ’° d2 


-3 


P + 5 pB ,EoT 
P +7 p 7 , o 


(28) 


(1 + z v d) 


o P + 5 
* P+7 


(1 + zf 


P+4 

P+7 


To calculate pb,e oT we consider the era when turbulence 
is generated by photons after neutrino decoupling. As 
already argued in the Section |II A[ we can trace the evo¬ 
lution of the magnetic energy all the way back to the time 
of magnetogenesis using the turbulent decay law given in 
Eq. §. This is possible since the turbulent stage due to 
neutrinos followed by the viscous stage before neutrino 
decoupling can be well approximated by only consider¬ 
ing a turbulent damping stage throughout the epoch T5]. 
Hence, using the turbulent damping decay law, we can 
write 


PB, EoT = Pb,* 


) ?£+Z 

P+7 


(29) 


t DC {z) 2.06 x 10 33 
“ 0+2 


(I-IYpY'z 


(24) 


In the cases where the viscous regime starts before the 
p-e ra, i.e. TeoT > 7+i, we can use Eqs. (281 and (291 to 
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estimate the p-type distortion from Eq. (25), we find 
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n P+5 


M = 


5 P + 7 V Py,o ) \ T* ) 


X 


rZend 2 P+ 4 

/ dz(l + 2 ) *>+ 7 exp 

«/ 2 ,; 


-(—) 
V^DC/ 


(30) 


where we integrate from Zi = z v d to z en d- In the cases 
where the turbulent regime has not ended by the start of 
the /i-era, e.g. applicable for fully helical magnetic fields 
generated at the QCD scale with £ ~ 1 — 10~ 2 as can be 
seen in figure |T] we find 


P = 


14 p + 5 / p B> * 


15 p + 7 VPt.o 


o P+5 
P+7 


(31) 


P+3 

dz (1 + z) p+7 exp 


\zncj 


where we integrate from Zj to z enc i = zeoT- 

Equations (30) and (31) above are valid for either non¬ 


helical fields (p = 0) or maximally-helical fields (p = —4). 
Since we are interested in different initial conditions, in 
particular varying initial helicity fractions /*, we must 
consider the full evolution history to calculate the spec¬ 
tral distortions. This is done numerically and the results 
are shown in Figures [2] and [3] for magnetic fields gener¬ 
ated at the EW and QCD scales respectively. The p- 
type spectral distortion varies with p = p (e, T*, /*). In 
order to maximise this distortion we can set the maxi¬ 
mal value e = pb,*/Pj,o ~ 1 at the time of magnetoge¬ 
nesis, which corresponds to an in itial field strength of 
B\* — 3 x 10” 6 G. From Eq. (30) we see that if e < 1 


then TeoT will be larger, thereby increasing the chemical 
potential. However, Te 0 t does not depend so strongly on 
e, therefore £ should be maximised in order to maximise 
p. Hence, with e = 1, integrating Eq. p0| we obtain 
upper limits from Figures [2] and [3] For fields generated 
at the EW phase transition T* ~ 100 GeV, see Figure [5] 
with /* < 10“ 14 we find 


|p| < 1 x 10" 10 . (32) 

The above satisfies the current COBE/FIRAS limit 
|p| < 9 x 10 -5 [2], but will also not be detectable by 
a new PIXIE-like experiment which would place a new 
upper limit of |p| < 5 x 10 -8 d?j. From this we con¬ 
clude that causally generated non-helical magnetic fields 
at the electroweak phase transition will not produce any 
detectable CMB p-type spectral distortions. This is true 
even if we consider the non-helical inverse transfer effect 
seen in Refs. [M1I4T] and discussed below Eq. ([5]), see the 
gray dashed line in Figure [2] For non-helical fields gen¬ 
erated at the QCD phase transition T* ~ 200 MeV, see 
Figure[3j we obtain the upper limit |p| < 3 x 10 -8 , which 
satisfies the current COBE/FIRAS limit [T4], and is very 
much on the limit of detectability by a new PIXIE-like 
experiment m The results do not change much if we 


Log 10 (e) 

-2.0 -1.5 -1.0 -0.5 0.0 



FIG. 2. EW scale: In plot (a) we show the p-type distortion 
generated due to the decay of magnetic energy initially gen¬ 
erated at the EW scale. Here we plot the spectral distortion 
p vs £, where £ = Pb,*/ P 7 ,o ~ 1 corresponds to an initial field 
strength — 3 x 10~ 6 G. The (solid, blue) lines from top 
to bottom, in both plots, correspond to initial helicity frac¬ 
tions = {10 —3 ,10 -4 ,10 -6 , < 1(T 14 }. The maximally he¬ 
lical case ft = 1 (solid, red) is also shown. In plot (b) we 
show the final field strength Bo and coherence length A b that 
would be observed today, i.e. after MHD turbulent decay, see 
Eq. \22\ . We also show the approximate constraint on mag¬ 
netic fields from CMB observations, Bq < 10 _9 G see Ref. [48j 
and references therein. The results for non-helical magnetic 
fields with an inverse transfer of energy (see discussion below 
Eq. §) are also shown (dashed, gray lines). 


consider the non-helical inverse transfer effect discussed 
below Eq. ([5]), see the gray dashed line in Figure [3] 

The situation is quite different if magnetic helicity is 
introduced. As can be seen from Figures [2] and [3} if the 
initial helicity fraction /* > (l0 -3 — 10 -4 ) then PIXIE- 
observable p-type distortions can be generated. We also 
note that for maximally-helical fields /* = 1 the current 
COBE/FIRAS limit is not broken )14j . so that we can¬ 
not constrain primordial helicity from current data. If a 
future PIXIE-like experiment positively detects a p-type 
spectral distortion then primordial magnetic helicity can 
be constrained. 

Here we note that in Refs. ESS ED authors considered 
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FIG. 3. QCD scale: In plot (a) we show the y-type distor¬ 
tion generated due to the decay of magnetic energy initially 
generated at the QCD scale. Here we plot the spectral distor¬ 
tion p vs e, where e = pB,*/p-y,o ~ 1 corresponds to an initial 
field strength B\, t — 3 x 10 -6 G. The (solid, blue) lines from 
top to bottom, in both plots, correspond to initial helicity 
fractions /* = {10 —1 , 10 -4 ,10 -6 , < 10~ 14 }. The maximally 
helical case /* = 1 (solid, red) is also shown. In plot (b) we 
show the final field strength Bo and coherence length A b that 
would be observed today, i.e. after MHD turbulent decay, see 
Eq. (22 l. We also show the approximate constraint on mag¬ 
netic fields from CMB observations, Bq < 10 _9 G see Ref. (IS. 
and references therein. The results for non-helical magnetic 
fields with an inverse transfer of energy (see discussion be¬ 
low Eq. ©) are also shown (dashed, gray lines). Here we 
comment on the seemingly strange behaviour of the plot for 
/* = 1CP 1 in the QCD case (also applicable to the non-helical 
inverse transfer case, gray dashed line). As £ decreases the on¬ 
set of the viscous freezing regime begins earlier. For a certain 
value of £ the start of the viscous freezing regime coincides 
with the start of the y-era, hence there is practically no mag¬ 
netic energy decay in this time and hence y becomes very 
small. As £ decrease further, the start of the viscous damping 
stage occurs within the y-era and a large p can once again be 
generated. 


the evolution of the photon diffusion scale, i.e. the Silk 
damping scale, as the evolution of the damping scale kd- 
Where we can write kj « kd with an equivalent evolu¬ 
tion. However in these works, the authors only con¬ 


sidered a scale-invariant Alfven velocity, corresponding 
to scale-invariant spectrum n = —3, this gives the much 
faster evolution of the integral scale fcj oc a -3 / 2 (for non- 
helical fields with p = 0) as seen in their paper [2D]. How¬ 
ever, the magnetic energy is still considered to evolve 
along the large scale spectrum, i.e. Bi oc k 5 / 2 , hence 
the authors find a much faster magnetic field decay rate 
Pb oc a -15 / 2 . With this decay rate, the chemical po¬ 
tential scales as p oc / dz(l + z) 13 / 2 e _ ( Z//zDC 4 5/ " as seen 
in their paper [2D]. This fast magnetic field decay rate 
leads to a huge overestimation of the generated chemical 
potential for the causally generated magnetic fields. In 
our study, as can be seen in Figures [2] and [3] the maxi¬ 
mum possible magnetic energy e = 1, which corresponds 
to B\,* — 3 x 10 -6 G at magnetogenesis, does not over¬ 
generate spectral distortions. However, in the work of 
Ref. [2D], an upper limit of ~ 10 -11 nG on the field 
strength of non-helical magnetic fields is obtained due 
to current spectral distortions constraints. 


B. CMB y-type distortions from decaying magnetic 

fields 


After the end of the p-e ra, for z < 5 x 10 4 , the elastic 
Compton scattering becomes ineffective. From then on 
the spectral distortion becomes a y-type i.e. defined by 
the Compton y-parameter ilB, 


_ 1 r dec dz dQ 
V 12 Jz end Pi dz ’ 


(33) 


where z eq — 3265 and z^ec — 1088 for the time of matter- 
radiation equality and decoupling respectively. It is pos¬ 
sible that intermediate distortions between p and y-type 
are generated [25) [28], but for the purpose of this paper 
we will only calculate the p and y defined above. 

We can see from the plots in Figure [l] that in all cases, 
for magnetic fields generated at the EW or QCD scale 
and for any initial helicity fraction, the plasma is in a 
viscous state (viscous freezing or viscous damping) in the 
y-era. Hence, as a first approximation for the maximum 
y-type distortion generated we can assume the viscous 
damping law given in Eq. ([7]) and a radiation dominated 
universe throughout the y-era. From Eq. (33) and the 
above considerations we find 


P + 5 / p B} * \ f Teqt \ 2 
4 P + 7 V Pi, o ) \ T * J 


x (1 + z v d) 3 "+ 7 


/■ 

J Z P 


* ^ o pm-* 

d£ (1 + z) p + 7 


(34) 


Again, the above expression is valid for non-helical fields 
(jp = 0) or maximally helical fields (p = — 4). For varying 
initial helicity fractions /* a full calculation is done taking 
into acound the full evolution history. The results are 
shown in the plots of Figure [4] 























For non-helical fields, we can see from Figure [4] 
that the maximum possible y-type distortion gener¬ 
ated, corresponding to fields generated at the QCD 
scale, is y < 8 x 10 -10 . This result satisfies the current 
COBE/FIRAS limit y < 1.5 x 10 -5 [H], and is probably 
not detectable by a new PIXIE-like experiment, which 
would place a new lower bound at y < 10 -8 if not de¬ 
tected El- For non-helical fields generated at the EW 
phase transition there is little hope for detection with 
the maximum distortion at y ~ 3 x 10~ 14 . We note that, 
if we consider the non-helical inverse transfer effect dis¬ 
cussed below Eq. |5]) [see the gray dashed line in Figure 
[|, the y-type distortion is still undetectable for fields 
generated at the EW scale, however it seems now possi¬ 
ble to detect a y-type distortion if the fields are generated 
at the QCD scale. 

For decaying helical magnetic fields there is a possi¬ 
bility to detect the y-type distortion. However, since the 
y-e ra is later in the evolution history, magnetic fields have 
already substantially decayed and thus generate a smaller 
signal than that for the /r-type distortion. Hence y-type 
distortions will always provide weaker constraints than 
/r-type distortions. 


IV. CONCLUSIONS 

Magnetic fields generated in the very early Universe 
decay in the radiation dominated epoch due to turbulent 
MHD effects. The decaying magnetic fields inject en¬ 
ergy into the primordial plasma which can lead to /i-type 
and y-type distortions to the CMB black body spectrum. 
The current COBE/FIRAS limits on these spectral dis¬ 
tortions are very tight |/i| < 9 x 10~ 5 and y < 1.5 x 10~ 5 
P~4). However there is the exciting possibility of a new 
PIXIE-like experiment which could place much stronger 
upper limits of |y| < 5 x 10~ 8 and y < 10 -8 if no detec¬ 
tion is made (4JJ . Any prediction for spectral distortions 
above the PIXIE limits is what we call detectable. 

In this paper we consider the evolution of helical and 
non-helical magnetic fields generated by some causal pro¬ 
cess in the early Universe. We calculate the spectral dis¬ 
tortions using the decays laws of Refs. [HH31I33]. We 
find that causally generated non-helical magnetic fields, 
with an initial helicity fraction less than ~ 10 —14 , gener¬ 
ated at the EW phase transition will not produce any de¬ 
tectable CMB y-type or y-type spectral distortions. This 
remains true even if the inverse transfer effect for non¬ 
helical fields seen in Refs. [31II] is considered. Hence, 
to produce observable spectral distortions from the decay 
of magnetic fields generated at the EW phase transition, 
a non-negligible helical component is required. 

Here we note that, if the inverse transfer effect for non¬ 
helical fields is applicable Si], it looks possible to gen¬ 
erate small amounts of detectable distortions from mag¬ 
netic fields generated at the QCD phase transition. We 
also note that magnetogenesis at the QCD phase tran¬ 
sition is disfavoured compared to magnetogenesis at the 


Log 10 (£) 

-2.0 -1.5 -1.0 -0.5 0.0 



FIG. 4. In plots (a) and (b) we show the y- type 

distortion produced due to the decay of magnetic en¬ 
ergy initially generated at the EW and QCD scales re¬ 
spectively. Here we plot the spectral distortion y vs 
e, where e = Pb,*/p~/,o ~ 1 corresponds to an initial field 
strength of B\,, — 3 x 10 -6 G. In plots (a) and (b) the 
(solid, blue) lines from top to bottom correspond to ini¬ 
tial helicity fractions /* = {lO -3 ,10 -4 ,10 -6 , < 10 -14 } and 
/* = {10 —1 ,10 —4 ,10 -6 , < 10~ 14 } respectively. The maxi¬ 
mally helical case f* = 1 (solid, red) is also shown. The final 
field strength Bo and coherence length A b that would be ob¬ 
served today, i.e. after MHD turbulent decay, see Eq. ( |22| , 
are the same as in plots (b) of Figures [5] and [3] We also 
show (dashed, gray lines) the results from non-helical mag¬ 
netic fields with an inverse transfer of energy (see discussion 
below Eq. ©)■ 


EW phase transition. Under early Universe conditions 
with very small chemical potentials the QCD phase tran¬ 
sition is a smooth transition [35] whereas the EW phase 
transition could be first-order in certain Standard Model 
extensions m- 

The conservation of magnetic helicity in the early Uni¬ 
verse leads to an inverse cascade of energy and the slow¬ 
ing down of magnetic decay for fully helical fields. This 
means that, at the time when CMB spectral distortions 
can be generated, the magnetic field amplitude is rel- 
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atively large compared to the non-helical case. This 
can lead to the generation of larger spectral distortions. 
If CMB spectral distortions are observed by some new 
PIXIE-like experiment, then it is likely that magnetic 
helicity plays an important role. However, there is a de¬ 
generacy in the parameter space, since different parame¬ 
ter sets can give the same spectral distortions signal. For 
example, fields generated at the QCD phase transition 
with smaller e = Pb,*/ jo 7 ,o and/or helicity can produce 
the same //-type distortions as fields generated at the 
EW phase transition but with larger e and/or helicity. 
However, if a /t-type distortion is detected by a PIXIE- 
like experiment, it would rule out non-helical magnetic 
fields produced at either the EW or QCD phase transi¬ 
tion. A positive detection would give us a lower bound 
on the primordial magnetic helicity. The lower bound 
would be somewhere of the order /* > (10 -4 — ICC 3 ). 
This is much greater than the primordial magnetic he¬ 
licity generated in the simplest models of EW baryoge- 
nesis [ST; where /* ~ 10 -24 assuming B \= -B™* x and 
As,* = Aew -523 • However, there are new mechanisms 
being proposed recently which can excite magnetic he¬ 
licity in the early Universe due to the Chiral anomaly 
[531 El. It will be interesting to investigate such mecha¬ 
nisms in the future. 


It is also interesting to mention the recent tentative 
observations of large scale helical magnetic fields from 
7 -ray observations [55]. Such studies have seen some evi¬ 
dence, albeit rather weak, of fully helical fields of strength 
10~ 14 G on scales of 10 Mpc. If such fields originated 
from a time before the /t-era, then it is possible that 
such observations would be accompanied by a detectable 
signal for a PIXIE-like experiment. The combination of 
such two observations would be compelling evidence for 
large scale helical magnetic fields. We also note that the 
CMB distortions anisotropies (see e.g. Refs. [221 155)i. 
albeit potentially very hard to detect, could give inter¬ 
esting signals due to the large helicity of the magnetic 
fields. Unique signatures in the spatial correlations are 
expected due to the helical nature of the magnetic fields. 
This will be investigated in future publications. 
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